Seismic radiation from dynamic coalescence, and the reconstruction of dynamic source parameters on a planar fault     

Nobuki Kame  Koji Uchida 《Geophysical Journal International》2008,174(2):696-706

The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω⁻²  -type displacement spectra.  相似文献   

Dynamics of continental collision: influence of the plate contact     

Roberta De Franco  Rob Govers  Rinus Wortel 《Geophysical Journal International》2008,174(3):1101-1120

  相似文献   

Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modelling   总被引：1，自引：0，他引：1  

Hongwei Gao  Jianfeng Zhang 《Geophysical Journal International》2008,174(3):1029-1036

The perfectly matched layer (PML) absorbing boundary condition is incorporated into an irregular-grid elastic-wave modelling scheme, thus resulting in an irregular-grid PML method. We develop the irregular-grid PML method using the local coordinate system based PML splitting equations and integral formulation of the PML equations. The irregular-grid PML method is implemented under a discretization of triangular grid cells, which has the ability to absorb incident waves in arbitrary directions. This allows the PML absorbing layer to be imposed along arbitrary geometrical boundaries. As a result, the computational domain can be constructed with smaller nodes, for instance, to represent the 2-D half-space by a semi-circle rather than a rectangle. By using a smooth artificial boundary, the irregular-grid PML method can also avoid the special treatments to the corners, which lead to complex computer implementations in the conventional PML method. We implement the irregular-grid PML method in both 2-D elastic isotropic and anisotropic media. The numerical simulations of a VTI lamb's problem, wave propagation in an isotropic elastic medium with curved surface and in a TTI medium demonstrate the good behaviour of the irregular-grid PML method.  相似文献   

Interpolation using a generalized Green's function for a spherical surface spline in tension   总被引：1，自引：0，他引：1  

P. Wessel  J. M. Becker 《Geophysical Journal International》2008,174(1):21-28

A variety of methods exist for interpolating Cartesian or spherical surface data onto an equidistant lattice in a procedure known as gridding. Methods based on Green's functions are particularly simple to implement. In such methods, the Green's function for the gridding operator is determined and the resulting gridding solution is composed of the superposition of contributions from each data constraint, weighted by the Green's function evaluated for all output–input point separations. The Green's function method allows for considerable flexibility, such as complete freedom in specifying where the solution will be evaluated (it does not have to be on a lattice) and the ability to include both surface heights and surface gradients as data constraints. Green's function solutions for Cartesian data in 1-, 2- and 3-D spaces are well known, as is the dilogarithm solution for minimum curvature spline on a spherical surface. Here, the spherical surface case is extended to include tension and the new generalized Green's function is derived. It is shown that the new function reduces to the dilogarithm solution in the limit of zero tension. Properties of the new function are examined and the new gridding method is implemented in Matlab® and demonstrated on three geophysical data sets.  相似文献   

三种猪毛菜种子耐盐性与幼苗干燥存活能力   总被引：3，自引：1，他引：2  

刘鹏  田长彦  冯固  薛英 《干旱区地理》2008,31(2):271-277

在控制条件下研究了三种藜科植物紫翅猪毛菜(Salsola affinisC.A.Mey.)、浆果猪毛菜(Salsola foliosa(L.)Schrad.)和钠猪毛菜(Salsola nitrariaPall.)种子的萌发特性以及萌发后的幼苗在室温中放置14天和30天之后重新放置在原溶液中的存活状况。结果表明:三种猪毛菜种子在无菌水、0.085 mol/L、0.17 mol/L、0.34 mol/L NaCl溶液中萌发的很好,但是在0.51 mol/LNaCl中萌发率很低,三种猪毛菜幼苗在干燥30天重新放置在原溶液中以后仍然具有很高的存活率。上述结果表明:已经萌发了的幼苗在干燥环境中仍然保持存活的能力是盐生植物为了适应极端干旱和盐渍荒漠环境的一个重要的生存策略。  相似文献   

一种基于模糊正域类的模糊属性约简方法     

冯楠坪  周磊 《成都信息工程学院学报》2013,28(4):434-439

为了实现模糊信息系统属性约简问题,探讨一种基于正域类的模糊属性约简方法。利用模糊二元相似关系建立模糊相似类的方法,进而定义决策属性关于条件属性的模糊正域类,得到正域类属性的重要度。结果表明:一方面,实现模糊信息系统属性约简;另一方面,通过取不同的相似精度,得到不同属性约简集。  相似文献   

Focal Mechanism Solutions of the 2008 Wenchuan earthquake sequence from P-wave polarities and SH/P amplitude ratios: new results and implications     

Yuan Tian  Jieyuan Ning  Chunquan Yu  Chen Cai  Kai Tao 《地震科学（英文版）》2013,26(6):357-372

The 2008 Wenchuan earthquake, a major intraplate earthquake with M _w 7.9, occurred on the slowly deforming Longmenshan fault. To better understand the causes of this devastating earthquake, we need knowledge of the regional stress field and the underlying geodynamic processes. Here, we determine focal mechanism solutions (FMSs) of the 2008 Wenchuan earthquake sequence (WES) using both P-wave first-motion polarity data and SH/P amplitude ratio (AR) data. As P-wave polarities are more reliable information, they are given priority over SH/P AR, the latter of which are used only when the former has loose constraint on the FMSs. We collect data from three categories: (1) permanent stations deployed by the China Earthquake Administration (CEA); (2) the Western Sichuan Passive Seismic Array (WSPSA) deployed by Institute of Geology, CEA; (3) global stations from Incorporated Research Institutions for Seismology. Finally, 129 events with magnitude over M _s 4.0 in the 2008 WES are identified to have well-constrained FMSs. Among them, 83 are well constrained by P-wave polarities only as shown by Cai et al. (Earthq Sci 24(1):115–125, 2011), and the rest of which are newly constrained by incorporating SH/P AR. Based on the spatial distribution and FMSs of the WES, we draw following conclusions: (1) the principle compressional directions of most FMSs of the WES are subhorizontal, generally in agreement with the conclusion given by Cai et al. (2011) but with a few modifications that the compressional directions are WNW–ESE around Wenchuan and ENE–WSW around Qingchuan, respectively. The subhorizontal compressional direction along the Longmenshan fault from SW to NE seems to have a left-lateral rotation, which agrees well with regional stress field inverted by former researchers (e.g., Xu et al., Acta Seismol Sin 30(5), 1987; Acta Geophys Sin 32(6), 1989; Cui et al., Seismol Geol 27(2):234–242, 2005); (2) the FMSs of the events not only reflected the regional stress state of the Longmenshan region, but also were obviously controlled by the faults to some extent, which was pointed out by Cai et al. (2011) and Yi et al. (Chin J Geophys 55(4):1213–1227, 2012); (3) while the 2008 Wenchuan earthquake and some of its strong aftershocks released most of the elastic energy accumulated on the Longmenshan fault, some other aftershocks seem to occur just for releasing the elastic energy promptly created by the 2008 Wenchuan earthquake and some of its strong aftershocks. (4) Our results further suggest that the Longmenshan fault from Wenchuan to Beichuan was nearly fully destroyed by the 2008 Wenchuan earthquake and accordingly propose that there is less probability for great earthquakes in the middle part of the Longmenshan fault in the near future, although there might be a barrier to the southwest of Wenchuan and it is needed to pay some attention on it in the near future.  相似文献   

2010年6月广西凌云-凤山交界3级震群应力水平分析     

郭培兰  龙政强  阎春恒  陆丽娟  牟剑英 《华南地震》2014,(1):11-20

2010年6月28日至8月15日,在广西凌云与凤山交界持续发生了高密集、高集中、高频率的震群活动,共记录地震3 071次,其中ML2~2.9级地震41次,3级以上地震3次,最大为7月1日10时27分3.2级地震。对此次震群的空间分布范围、2级以上地震震源机制解P轴优势方向、b值特征、应力降等进行了分析研究。结果表明,此次震群震中集中,深度较浅,震源机制解P轴没有形成优势分布,序列参数b值较高,约为1.21,应力降主要分布在0.19~2.0 Mpa之间,认为广西凌云-凤山交界3级震群是属于低应力背景下的地震活动,不具有前兆震群和前震序列的性质。分析结果对正确了解此次震群的特征、正确判断地震趋势具有一定的帮助作用。  相似文献   

Analytical model of well leakage pressure perturbations in a closed aquifer system     

《Advances in water resources》2014

Deep saline aquifers are commonly used for disposal and storage of various surface fluids. The target injection zone must be hydraulically isolated from overlying zones in order to ensure containment of the injected fluids. Improperly plugged nonoperational abandoned wells that penetrate the injection zone are the main potential leakage pathways. Leakage through such wells may cause an observable pressure signal in a zone overlying the injection zone; such a signal can be used to detect the leakage. In this paper we develop an analytical model to evaluate the pressure change induced by leakage through a well in a multilayer system. Unlike previous analytical models on the topic, our model uses a closed system, which may significantly affect the strength and behavior of the pressure signal induced by leakage. The analytical model is first presented for a two-layer system centered at the leaky well location. We evaluate the leakage-induced pressure change using the Laplace transform of Duhamel’s superposition integral, yielding the solution in the Laplace domain. We then derive a late-time asymptotic solution using the final value theorem, which suggests that the leakage rate becomes constant after sufficient time. We then obtain the multilayer solution by extending the two-layer solution and presenting it in matrix form in the Laplace domain. We apply the solution to three examples. In the first example, we apply the analytical model to a two-layer system, investigating its behavior and comparing the results with a numerical solution. In order to demonstrate behavior and potential applications of the multilayer analytical model, we present two multilayer examples: one with identical layers and another, replicating a CO₂ storage site, with dissimilar layers. The leakage-induced pressure change does not necessarily decrease as the distance increases from the injection zone toward the surface.  相似文献   

Comparison of alternative upscaled model formulations for simulating DNAPL source dissolution and biodecay     

Jack C. Parker  Ronald W. Falta 《Advances in water resources》2008

This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution.  相似文献